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SUMMARY 


An unsteady finite difference procedure is used to calculate the steady 
inviscid flow over a 11^6 thickness ratio NASA supercritical profile of 

LWP 505 at M =0.80 and cr ■ 0. An attempt is made to include the viscous 

00 

effects using a modified form of Head's entrainment method to calculate the 
turbulent boundary layer. The results of the inviscid calculations compared 
to the experimental results of Whitcomb indicate a significant influence of 
viscosity, though the comparison is partially obscured by the unknown wind 
tunnel wall interference effects. The effects of viscosity were manifested 
in three ways. Firstly, the experiments clearly indicated the presence of 
shock-induced separation at the base of the shock. The resultant flow dis- 
placement effect produced a significant upstream movement of the shock. 
Secondly, the displacement effect of the boundary layer, primarily on the aft 
upper surface, reduced the effective eft camber which in turn resulted in 
changes of the "plateau" pressures on both the upper and lower surfaces charac- 
teristic of a decreased flap effect. Thirdly, the severe adverse gradients 
on the aft portion of the lower surface resulted in a significant thickening 
or a local separation of the boundary layer that led to decreased over- 
pressures. 

The attempt to predict the viscous effects using the compressible form 
of Head's integral method with a modified auxiliary equation f the form 
factors was unsatisfactory. Though a reasonable separation buuule was es- 
tablished on the lower surface, a grossly exaggerated displacement effect 
resulted downstream of the shock on the upper surface. Here in these explora- 
tory calculations the experimental pressure distribution was simply used to 
check the boundary layer equations. There clearly is substantial further 
effort required to evolve a satisfactory boundary layer procedure, which must 
then be coupled in a still unproven manner with the inviscid procedure. 
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INTRODUCTION 


The calculation of the planar supercritical flow over airfoils is uni- 
versally recognized as a horrendous task. Sources of the difficulties are 
manifold. At the high Reynolds numbers of pertinence the resulting boundary 
layer is turbulent, and is sufficiently thin that the classical boundary 
layer concept would be applicaole were it not for the presence of the termi- 
nating shock wave. The presence of the shock wave can 30 drastically alter 
the ’'effective" shape of the profile through the displacement effects of the 
boundary layer that the inviscid solution no longer serves directly as a 
viable basis to obtain the pressure distribution. In short the inviscid arid 
boundary layer flows are strongly coupled. Despite the strong coupling with 
proper ingenuity it would not be a hopeless task to compute the coupled flow 
if methods were available to calculate both the inviscid flow over a given 
shape, and the turbulent boundary layer flow, possibly separated, in the 
presence of a prescribed pressure distribution. 

In the present study the objective will be to carry out the inviscid 
calculations for a prescribed airfoil, namely the 11^6 thickness ratio NASA 
supercritical profile of LWP 505, at or = 0 and M = 0.80, and attempt to 

OD 

extend the results to include the viscous effects. 


THE INVISCID CALCULATIONS 


To compute the inviscid flow the vu.iteady finite difference procedure 
as described in Ref. 1 will be used. Here the coexistence of subsonic and 
supersonic flows with their highly differing mathematical characters in the 
steady state makes the unsteady approach more tractable than the steady from 
a calculational point of view, since in the latter representation the character 
of the equations remains hyperbolic in both the subsonic and supersonic flow 
regimes. 

In Figure 1 is shown the resulting Mach wave pattern for the supercritical 

profile at M = 0.80 and a = 0 , whereas in Figure 2 the resulting chordwise 
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distribution of the pressure is shown. Also shown in Fig. 2 are some experi- 
mental results as obtained by Whitcomb (Ref. 2) for several geometric values 
of a. The pressure distributions for the several values of a are shown 
here since the precise value of the angle of attack is not known because of 
the disturbing influence of the wind tunnel walls. If tests had been carried 
out additionally for a geometric a of approximately 2 , the plateau values 
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of the pressure on both the upper and lower surfaces would agree closely with 

the calculated result. The Immediate conclusion might then be that there is 

o 

a minus 2 correction in or due to wall effects for this case. That this 
may not be a completely correct conclusion may be suggested by the lowness 
of the trailing edge pressure, and the fact the shock pressure rise falls 
significantly short of that for the normal shock. The latter clearly indicates 
the presence of a separation bubble immediately downstream of the shock which 
then leads to a significant boundary layer displacement effect in the vicinity 
of the upper surface trailing edge as manifested by the depressed trailing 
edge pressure. The added displacement effect of the boundary layer here gives 
rise to a lessening of the effective trailing edge angle; that is, a lessening 
of the aft camber. This in turn leads to increased plateau pressures on the 
upper side and decreased plateau pressures on the lower side as would result 
from a decreased flap effect. In short, the modifications of the pressure 
distribution due to viscous effects on the airfoil surface are of the same 
general character as that due to wall interference effects. 

Aside from the viscous effects just described a more serious consequence 
of the separation bubble is that the shock wave at the surface is now oblique 
instead of normal as in the unseparated case, and the resultant shock pressure 
rise will then be significantly less than the normal shock pressure rise. As 
a result the shock wave will move upstream to a position sufficiently fa r up- 
stream to allow a sufficient "run" of su’ sonic recompression to attain the 
required trailing edge pressure. Quite apart from the above effect there is 
a countering influence on the shock movement by the depression of the trailing 
edge pressure. Thus, for a depressed trailing edge pressure there would be 
less subsonic recompression required, and as a result the shock wave would not 
have to migrate as far upstream as for a higher trailing edge pressure. 

Finally in Fig. 2 it is seen that there is a significant depression of 
the overpressures on the lower aft surface caused by a separation or a thicken- 
ing of the boundary layer resulting from the severe adverse subcritical pressure 
gradient caused by the aft camber. 

The above results clearly indicate that inviscid results must be supple- 
mented by the addition of the viscous effects. 


VISCOUS CALCULATIONS 


There presently does not exist a suitable procedure to calculate a 
turbulent boundary layer in the presence of a prescribed external pressure 
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distribution containing adverse pressure gradients sufficient to lead to 
separation. An attempt will be made however to extend an existing method for 
this purpose. Of the many procedures one of the more promising methods is 
due to Head (Ref. 3) and Green (Ref. 4) , and this procedure will be adopted 
for our further consideration. Another widely used procedure is due to Nash 
and Macdonald (Ref. 5)* The latter method is an extremely simple procedure 
and is based upon the use of the classical streamwise momentum integral 
equation which results in the following ordinary differential equation for 
the momentum thickness 0: 


dx 


- (H + 2 - M ) 0 
e 


d In u 
dx 


- + - C 

2 f 


( 1 ) 


Here H is the usual form factor H = 6 /f) , M and u are the Mach number 

6 6 

and velocity in the external inviscid flow, and C is the wall skin friction 

coefficient. In the Nash-Macdonald method Eq. (l) is supplemented by two 

additional equations in order to obtain a fully determined system, assuming 

for the time being that u (x) and M (x) are prescribed. First of the 

0 0 

required equations expresses as a function of the Reynolds number based 

upon 0, and one or more of the other dependent variables appearing in (l). 

A suitable form is that given in Ref. 5> or the well known Ludwieg-Tillmann 
formula in one of its many variant forms. The second auxiliary relation 
suggested by Nash and Macdonald is based upon the assumption that the flow 
is locally an "equilibrium boundary layer;" that is one that fulfills the 
condition 


G(n) 


6.1 (rr + 1.8l) 


1/2 
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1 + 0.178 M 


- 1 


and 
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tor a given u (x) and M (x) the differential equation (1) for 0 can 
0 . 0 

be integrated using ( 2 ) and the skin friction law, and one can obtain quite 
simply the chordwise distribution of 0(x), C f (x), as well as fl (x). For 
an unseparated boundary layer use of the above procedure has been found to 
yield a s ati s f actory distribution for 0 (x), but is inadequate to yield an 
acceptable distribution for 6 "(x) (Ref. 5). (The general experience among 
most integral methods has been that all procedures appear to yield a reasonable 
distribution for 0 , but the determination of ft* is a much more delicate 
undertaking.) For the separated case it can be expected that the Nash- 
Macdonald procedure would simply be grossly inadequate to predict fi (x). 

(Note that the authors of Ref. 5 do not claim the applicability of their 
procedure to the separated case) . 

A more sophisticated approach is the use of the compressible form of 
Head's procedure (Ref. 3) given by Green (Ref. 4). The basic equations in 
this procedure in addition to Eq. (l) include the integral form of the conti- 
nuity equation, namely, 


dA 

dx 


? 

F + (M 
e e 


1) A 


d In u 
e 

dx 


(3) 


Here 

AM • «' 


F = - — + ^ (Entrainment function) 

e u dx 

e 

where p is the density, 6 the boundary layer thickness, v the transverse 
velocity, and the subscript e denotes the local inviscid value at y = fi. 

In the Head-Green method the basic differential equations (l) and ( 3 ) 
must be supplemented by three additional equations. These are first a skin 
friction law, for which we use the Felsch modification of the Ludwieg-Tillmann 
law. (Ref. l). The second is an expression for the entrainment function F , 
and the third is a relation of the form H = H^) where H 5 a/ 0. It is 6 
in the last relation that an attempt was made to derive an expression valid 
for the separated case using available experimental data. For this purpose 
the experiments of Green (Ref. 4) and Seddon (Ref. 6 ) were used to derive the 


u 
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dy , and 
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relation H ■ H(H ) despite the fact that in the latter experiments it can 
be clearly shown that there is serious departure from the planar condition 
due to the influence of the sidewall boundary layer displacement effects 
(see Ref. l). 

To test the resulting set of equations calculations were carried out where 

an existing experimental pressure distribution from LWP 505 was used. The 

results showed a reasonable fi* upstream of the abrupt adverse gradients 

on both the upper and lower surfaces as would be expected. On the lower 

surface the abrupt recompression led to a local separation bubble that appeared 

to be reasonable, but downstream of the shock on the upper surface the resulting 
* 

distribution of 6 grew excessively towards the trailing edge. 

This unreasonable growth of 6 * as the trailing edge was approached 
could be directly traced to the inadequacy of the auxiliary relation H = HfH^) 
used which yielded near-singular values of H for the values of that 

were being predicted. 

Despite this shortcoming the resulting 6 *, with the unreasonable values 
arbitrarily modified, was used to obtain a modified airfoil shape, and the 
inviscid flow was recalculated in order to obtain a qualitative effect of 
the boundary layer displacement. The results indicated the change of the 
pressure distribution described earlier for the decreased aft camber. 

These highly tentative attempts to incorporate ths viscous effects clearly 
point out areas of future effort required to improve the procedure. With a 
reasonable amount of additional effort it is not excessively optimistic to 
expect that a viable boundary layer procedure can be eventually evolved to 
handle the specialized case of shock-induced separations on supercritical 
airfoil flows. 


Convair Division of General Dynamics 

San Diego, California, 18 September 1970 
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Figure 1. Mach Wave Pattern for the Basic NASA Profile 
at = 0. 80 and a = 0. 0° 
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